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INTRODUCTION 


Block,  von  Bloh,  and  Schellnhuber  (BBS)  (ref  1)  presented  a  box-counting  algorithm  fcv  the 
measurement  of  generalized  (fractal)  dimensions,  which  they  refer  to  as  "Efficient  Box-Counting"  (EBC). 

Meisel,  Johnson,  and  Cote  (ref  2)  described  a  box-counting  algorithm,  which  is  particularly  suited  to  the  analysis 
of  experimental  data  with  inherent  uncertainty  and  large  data  sets,  which  they  referred  to  as  "Agglomeration 
Box-Counting"  (ABC).  Atmanspacher,  Scheingraber,  and  Wiedenmann  (ref  3)  applied  techniques  similar  in 
spirit  to  EBC  and  ABC,  which  they  referred  to  as  "histogram  techniques,"  and  a  technique  based  on  the 
generalized  correlation  integral  (ref  4),  to  the  determination  of  the  muldhactal  properties  of  the  distribution  of 
galaxies  on  the  celestial  sphere  (ref  3). 

It  was  shown  in  Reference  2  that  although  box-counting  algorithms  converge  to  analytic  results  for 
model  fractal  point  sets  at  q  >  0,  they  fail  to  converge  to  analytic  results  at  q  <  0.  Similarly,  Atmanspacher 
et  aL  (ref  3)  reported  consistent  results  for  the  left  branch  of  the  f(a)  spectrum  (q  >  0),  but  noted  discrepancies 
for  the  right-hand  branch  (q  <  0).  It  is  evident  that  a  reliable  algorithm  for  the  determination  of  D(q)  at  q  <  0  or 
equivalently  for  the  determination  of  the  right-hand  branch  of  the  f(a)  spectrum  is  urgendy  required. 

The  present  work  describes  a  Box-Based  Correlation  Integral  (BBCI)  method  for  multihactal  analysis. 
The  BBCI  algorithm  has  been  tested  for  a  selection  of  model  multifractal  point  sets,  which  are  subject  to  simple 
analysis  along  the  lines  in  Halsey  et  al.  (ref  S).  BBCI  values  of  D(q)  converge  close  to  analytic  values  at  all 
q  e  [-25,25]  for  the  model  multifiactal  point  sets.  The  algorithm  is  well-suited  to  experimental  data  application, 
and  it  converges  at  least  as  rapidly  as  ABC  at  q  >  1. 

THEORY 

The  Box-Based  Correlation  Integral  (BBCI)  algorithm  is  implemented  as  follows: 

1.  Define  a  set  of  nested  hypercubes  ("boxes")  appropriate  for  the  given  (or  anticipated)  point  set  S. 

We  refer  to  these  boxes  as  elementary  hypercubes  or  elementary  boxes.  For  experimental  data,  the  appropriate 
choice  of  elementary  box  edge  should  reflect  the  inherent  uncertainty  of  the  coordinates  of  the  point  set 
within  limits  set  by  storage  requirements. 

Storage  requirements  are  effectively  determined  by  the  choice  of  elementary  hypercubes  rather  than  the 
size  of  the  point  set.  Thus,  large  point  sets  can  be  accommodated  with  relatively  modest  resources. 

2.  Compute  or  measure  the  occupation  numbers  nj  for  each  of  the  elementary  hypercubes.  Define  N  as 


ei  sawn  Cary  boxea 
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In  our  tests  N  ranged  between  about  1(H  and  10*. 

3.  Define  a  reference  set  that  comprises  a  subset  of  the  occupied  elementary  boxes.  Let  N,^  be  the 
number  of  elementary  boxes  in  the  reference  set.  In  our  tests  with  N  =  768*  elementary  boxes,  we  took  n„,  as 
the  smaller  of  the  to^  number  of  occupied  boxes  or  15,(X)0. 

4.  For  each  member  of  the  reference  set,  defme  a  set  of  hypercube  edge  lengths,  E  =  (2n+l)  Eo,  where 
n=0,l,2,...nvux-  In  our  tests  with  N  =  768*  elementary  boxes,  we  took  n^,^  =  24. 
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S.  For  each  q  of  interest: 

a.  Compute  the  box-based  generalized  correlation  integrals  C(qf  which  are  defined  as 

alaamaeazy  boxma  i 

C(g.F,£'o)  =  [i(  ^  n,njG{E,Eo,x^-Xj) 

for  the  E  values  defined  in  stq>  4,  where  r  runs  over  the  reference  set,  j  runs  over  aU  elementary  boxes, 

Urmt 

^2  =  E 


and  the  function 


G{E.Eq,x)  =  1,  if  \x^  \  {E-E^)/2,  for  all  components  of  x 

-  0,  otherwise, 

selects  the  elementary  boxes  contained  in  larger  hypeicubes  of  edge  E.  The  vectors  x,  and  Xj  point  to  elementary 
hypercubes  rather  than  members  of  the  point  set  Reducing  the  sum  over  elementary  boxes  in  the  side  E 
hypercube  by  unity  eliminates  self-correlation  contributions. 

The  box-based  generalized  correlation  integral  reduces  to  the  standard  generalized  correlation  integral 
(ref  4)  in  the  case  when  the  coordinates  of  the  members  of  the  multifiactal  subset  in  question  are  precisely 
known,  the  elementary  hypercubes  are  small  enough  to  contain  at  most  one  element  of  the  uiultifr^tal  subset, 
and  the  covering  set  is  comprised  of  hypercubes.  That  is, 

C(<J,E)  »  liinC(<3r/^,£'o)  • 
ab-o* 


Computation  was  slowed  down  and  convergence  was  not  improved  by  employing  (approximate) 
hyperspheres  rather  than  hypercubes.  Also  note  that  special  techniques  are  required  near  q  =  1;  here,  we 
restricted  q  such  that  iq  -  II  >  Vi. 

b.  Obtain  D(q)  and  the  corresponding  rms  error  in  D(q)  by  linear  regression  on  ln[C(qJ£,Eo)]  =  const  + 
D(q)  h^E]  for  E  e  ((2ih-1)  EqI  nsl,2,...nMAx}- 

RESULTS  AND  CONCLUSIONS 

The  BBCI  algorithm  has  been  tested  for  Euclidean  curves,  Koch  synunetric  and  asynunetric  triadic 
snowflakes,  split  snowflake  halls,  and  the  multifiactal  construction  based  on  Mandelbrot’s  13-element  generator. 
The  D(q)  determined  by  BBCI  converged  to  values  near  the  analytic  D(q)  values  for  q  tanging  b^een  -25  and 
■f2S  in  all  cases  as  N  was  increased. 
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Figure  1  shows  results  of  application  of  the  BBCI  algorithm  to  split  snowflake  halls  at  five  values  of  q. 
The  algcnithm  was  applied  to  a  randomly-oriented  construction  centered  on  768  by  768  square  boxes;  was 
the  lesser  of  15,000  or  the  total  number  of  occupied  boxes.  The  horizontal  lines  are  the  analytic  values. 
Convergence  was  within  about  2  percent  for  q  <  I  and  within  I  percent  for  q  >  1  in  this  test  Note  that 
convergence  is  about  an  order  of  magnitude  faster  for  q  >  1  than  it  is  for  q  <  1.  Except  that  convergence  is 
faster  for  constructions  with  lower  D(0),  Figure  1  is  typical  of  convergence  obtained  by  application  of  the  BBCI 
algorithm  to  the  constructions  tested. 

Figure  2  shows  the  results  of  application  of  the  BBCI  algorithm  as  in  Figure  1  for  N  »  2.14  10*  and 
analytic  values  of  D(q)  versus  q  for  split  snowflake  halls.  The  increase  in  the  error  in  D(q)  with  Iql  and  the 
order  of  the  discrepancies  are  apparent 

The  rms  errors  found  in  the  linear  regressions  were  of  the  order  of  the  differences  between  analytic  and 
BBCI  results  for  D(q).  An  expanded  discussion  of  the  results  of  application  of  BBCI  to  model  fractals  and  of 
the  effects  of  variations  of  number  of  elementary  boxes,  and  n^Ax  will  be  presented  in  a  fumre  publication. 
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Figure  1.  Analytic  and  BBCI  measurements  of  the  fractal  dimenuoa  D(q)  versus 
level  fOT  iandomlyK)riented  split  snowflake  halls.  (The  number  of 
points  in  the  construction  N  »  ll**^.) 
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Rgure  2.  Analytic  and  BBCI  values  of  the  fractal  dimension  D(q)  versus  q 
for  level  8  split  snowflake  hails.  The  line  is  the  analytic  curve. 
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